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ABSTRACT 

Current cc-SN models predict overproduction of 44 Ti compared to observations. We 
present a model for an alternative channel where a cc-SN explosion is followed by 
a neutron star detonation (Quark Nova or QN), resulting in a spallation reaction 
of SN ejecta that produces 44 Ti. We can achieve a 44 Ti production of ~ 1CP 4 M 
with our model under the right time delay between the QN and the SN. Our model 
also produces unique signals not found in standard, cc-SN nucleosynthesis models. 
Some of these unique signals include a significantly large production of 7 Be and 22 Na. 
We discuss some of these signals by analyzing the late time light curve and gamma 
spectroscopy of our model. 

Key words: spallation, core collapse supernovae, cc-SNe, Quark Nova, 44 Ti, nucle- 
osynthesis, SN1987A. 



1 INTRODUCTION 

Late time light curves of SNe are considered a relatively 
accessible resource for studying the nature of SNe, yet, the 
study of these light curves is ridden with challenges. Very few 
objects have been studied beyond ~ 200 days, because as the 
SN dims out, the dust and the background noise, makes the 
signals unreadable ( jLeibundgut fe Suntzcff 2003). Even the 
most extensively documented supernova, SN1987A, is dim- 
cult to study at its later stages, because it is very hard to 
extrapolate directly the bolometric luminosity after ~ 1400 



days (Leibundgut & Suntzeff 20031. The lack of late time 



light curve data make it difficult to formulate standard SN 
models. SN models have their parameters adjusted so that 
the isotopic yields reflect the abundances inferred from late 



time light curves (Limongi & Chieffi 2006), which makes 
them dependent on the quality and understanding of our 
current observations. An example of the problematic nature 
of standard core collapse SN models, is 44 Ti production. Al- 
though 44 Ti production is assumed to be universal in core 
collapse SN models ( |Woosley et"aL|l995[ ), only SN1987A is 
documented well enough to be able to extrapolate 44 Ti from 
the bolometric light curve. The problem is precipitated by 
the fact that current 7-ray telescopes hav e only been abl e 
to detect 44 Ti decay lines in Cassiopeia A ( The et al.|2006 1, 
whereas current models predict more detections. 
In this letter, we attempt to solve some of these problems 
by proposing an alternative channel to the standard cc-SN 
model. After the SN explosion, neutron stars that undergo 
an evolution (e.g. matter fallback, or spin-down) evolution 



towards quark deconfinement densities will transition to a 
lower energy, strange matter phase and experience an ener- 
getic detonation. (Quark Nova or QN) ( Ouyed et al.|2002 |. 
The neutron star will release its outer layers as relativistic, 
neutron rich ejecta, where the ejecta interacts with the SN 
envelope (double-shock Quark Nova, or dsQN). The time de- 
lay (hereafter idelay) between the two detonations constrains 
the nature of the interaction, where a time delay on the or- 
der of days will cause significant spallation of SN ejecta by 
QN neutrons. The spallation reaction creates unique signals 
in the late time light curve due to the modification of ra- 
dioactive isotope abundances. The time delay tdeiay between 
the QN and SN detonations acts as a parameter for the iso- 
topic abundances. The different abundances of isotopes will 
also create specific gamma spectroscopic signals. Currently, 
there has been many indirect observations for QN suggested 
(e.g. |Hwang fc Laming| | |20"l2) ; [Ouyed et al.| ( [20121 ); |Leahy| 



fc Ouyed (20081) 



In continuation to Ouyed et al. (20111, we explain the Ti 
synthesis measured in the late-time light curves of core col- 
lapse SN ( |Woosley fc Hoffman|1991 1 as a product of nuclear 
spallation from the dsQN, in contrast to the traditional view 
that it is produced by the SN in situ. The conditions nec- 
essary for the production of a QN, and the constraints set 
by the parameter tdciay, explains the rarity of observations 
of 44 Ti nuclear decay lines. In contrast to our earlier pa- 



per (Ouyed et al. 20111, we avoid excessive destruction of 
56 Ni by spallation, by taking into account the mixing of the 
SN's layers. Numerical studies have shown that Rayleigh 
Taylor instabilities can cause considerable mixing in a very 
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short time span (|Kifonidis et al.||2000{ |Hachisu et aL]p9T 



Mueller et al. 19911. Furthermore, observations in 1987A 



seem to su pport the picture of thorough mixing (Mueller 
et al.|l99l| . If 56 Ni is mixed through the SN envelope, there 



are less Ni atoms in the inner layer. This would lead to QN 
neutrons hitting less 56 Ni, as opposed to the picture where 
all 56 Ni is concentrated in the inner layers (Ouyed et al 



2011). By avoiding excessive nickel depletion, we can avoid 



the sub-luminosity characteristic in Ouyed et al. (20111. 



iVcoll = 



AR 
A 



%i,0.1 Mq ^ -^0, 0.1 Mq M C,Q.l Mq 



57.54- 



560-3 



16°. 3 



120-3 



^SAf.lOOOkm/s X t delay. 10 days 



(3) 

where A co n represents the number of collisions an incoming 
neutron with AR. The subscript indicates the units, for ex- 
ample, A/o,o.i m is in units of 0.1 Mq. The interaction of 
that neutron with a SN atom will release a multiplicity £ of 
nucleons, 



2 MODEL 



We extend the model we presented in |Ouyed et al.| ( |2011[ ) by 
including mixing. In analogy to lab terminology, we divide 
our model into a beam part, and a target part. The scenario 
consists of a beam of relativistic QN ejecta that collides 
with a target that is the expanding, inner layers of the SN 
envelope. The time delay between the SN denotation and 
neutron star detonation is the key parameter in our model. 



2.1 Beam 



According to past studies ( Keranen et al.|2005 Ouyed et al. 
20051, a detonating neutron star can eject its outer most 
layers as a mass of Mqn ~ 10 -3 Mq. We can model the QN 
ejecta as a pulse of A ~ 1.2 x 10 54 Mqn, -3 neutrons where 
Mqn, -3 means Mqn in units of 10 -3 Mq. The neutrons are 
relativistic, with an energy of ~ 10 GeV ( jOuyed et al.|2011| ). 



2.2 Target 

We model the target as an expanding SN envelope of inner 
radius R — i>sNtdeiay, where dsn = 5000 km/s is the velocity 
of the SN ejecta, and £deiay is the time delay between the neu- 
tron star detonation and the SN detonation. We assume the 
exploding star is ~ 20 Mq. In contrast to the original model 
( Ouyed et al.|20iT l, we assume mixing. The target is the in- 
ner Mi nn er ~ 1 Mq mass of the expanding SN envelope, 
which is composed of mixed 16 0, 12 C, and 56 Ni. We assume 
that whole SN envelope has a total 56 Ni mass of ~ 0.1 Mq, 
and O and C have total masses of ~ 1.0 Mq each ( |Ouyed| 



2012 ). If the mixing has spread each element homogeneously 



through whole SN envelope, then we can derive the masses of 
each element for the interior ~ 1.0 Mq target. The masses 
of the isotopes in the target are Mni = 0.0476 Mq, and 
Mo = M c = 0.476 M Q . We treat the target, total number 
density n tota i as constant, where, 



"•Total = n- N i + n Q + nc 



47ri? 2 A7?' 



(1) 



where AR is the thickness containing Mi nner and A to t a i 
is total number of atoms within AR. Our mean free path A 
is, 



A = 



1 



TlNiO-56 + Wc<Tl2 + WOO"16 ' 



(2) 



where we use the semi-emp irical cross section a a = 
45 mb x A 0,7 ( |Letaw et al |l983| ), where A is mass number. 
For our model, we use |1| and |2| to divide the thickness 
AT? into iV C oii imaginary layers of width A, 



((E, A) w 4.67^56(1 + 0.381n£)y n] 



(4) 



where 



20111 



1.25 < Y np < 1.67 (Cugnon et al. 1997 Ouyed et al 



2.3 Spallation Statistics 

At each imaginary layer k of radial thickness A, an atom will 
be bombarded by N^ its neutrons, resulting in the nucleus 
( Ouyed et aTpOn) , 



A 1 



A - t°(E°,A> 

3=0 



(5) 



Each element in the layer has a statistical weight Wi ~ 
Ni/ Ntotai ~ Mi/ (Ntotai AiTUH), where N is the total num- 
ber of atoms of isotope i, Mi is the inner mass of a specific 
isotope i, m,H is the mass of a proton, and Ai is the atomic 
number of isotope i. The statistical weight signifies the frac- 
tion of atoms of that particular element that will be hit by 
the incoming neutrons. 

TV^its and £o are drawn from a Poisson distribution that 
peaks at, 



A 



Anioni + Acctc + Nqoo 
N m + A c + A 



AT^l-e- 1 ) 
AttR 2 



(6) 

where Ai is the number of atoms of isotope i, and Oi is 
the spallation cross section for isotope i. 
For subsequent layers, it follows that. 



AL = (1 - e-V^Mts 



(7) 



and the mean nucleon energy for each subsequent layer 



is, 



E 



(8) 



where spallation ceases at an energy of E ~ 73 MeV 
( |Cugnon et al.|1997[ ). 

2.4 Light Curve Model 

We compute the late time light curve by extending our spal- 
lation code with the Cococubed program (Timmes). Our 
spallation code computes the mass yields of different iso- 
topes using the model in Section [2] and the Cococubed ex- 
tension computes the bolometric luminosity from the isotope 
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mass yields. We also take into account the contributions of 
57 Ni, and 60 Co to the bolometric light curve, and we take 



their mass values from Woosley et al. ( 1995). In our model, 



we assume the simple case where the bolometric luminos- 
ity is entirely caused by the thermalization of nuclear decay 
radiation. 



3 RESULTS 
3.1 Photometry 

The key isotopes for the light curve (Fig. [l} are 56 Ni and 
44 Ti. For tdciay = 4 days, our model synthesized ~ 10 -4 Mq 
of 44 Ti, which compares favourably to amounts found by 
7-ray telescopes in Cassiopeia A (Iyudin et al. 19941, 
and amounts estimated from the light curve of SN1987A 
(Suntzeff et~aI1|1992| [Woosley fc Hoffman||1991[). Further- 



more, we notice there is very little nickel that is depleted - 
indeed most of the original inner 0.0476 Mq is maintained 
for 3 days t delay ^ 10 days. If very little 56 Ni is depleted 
from the inner mass, then the total ~ 0.1 M© for 56 Ni that 
is spread throughout the whole envelope, is more or less 
maintained throughout different time delays. The ~ 0.1 Mq 
mass for 56 Ni compares quite well with what is observed 
in type II SN light curves ( |Suntzeff et al.||1992| |Woosley fc 



m ty 
|Hoffr 



Hoffman||l99l[ ). Another interesting isotope is Na, and it 
is synthesized at a high amount of almost ~ 1 0~ 4 Mq whic h 



contrasts with the estimates of ~ 10 Mq (Woosley et al 



19951 



This creates a higher contribution to the bolomet- 
ric luminosity, than normally expected. We see in the upper 
panel of Fig [l] a dependance of the light curve to tdciay • We 
can see that latter stages of the late time light curve are the 
most visibly affected by the time delay. This strong depen- 
dance in the later days is modulated by the production of 
44 Ti, because 44 Ti yields are very affected by the time delay 
(Fig.[3|. The earlier days of the late time light curve seem to 
be affected only weakly by the time delay change, because 
56 Ni destruction stays more or less constant (Fig. [3]). 



3.2 Isotope Abundance 

Our spallation model creates isotopes of almost every mass 
number in the range of A < 56. However, very few of the iso- 
topes created can be detected with 7-ray telescopes, because 
the SN envelope only becomes transparent to nuclear decay 



radiation after a time period of at least months ( Diehl 2012 1 , 
so only isotopes that have a half life of at least in the order 
of weeks will have detectable nuclear decay lines. Therefore, 
in Fig [3] we only include the masses of isotopes with half 
lives of at least in the order of weeks. 

In Fig. [2] and [3] we plot some of the trends of these 
isotopes. We notice that in Fig. [2j the plot heavily leans 
on light isotopes. This is caused by the spallation of mixed 
C and O in the inner layer. We also notice that the plot 
m Fig. [2] is divided in roughly two plateaus, where their 
boundary is in A ~ 16, which is the atomic mass number 
of Oxygen. In Fig. [3] we plot the mass of one of the lighter 
isotopes, 7 Be, which has a half life of 76 days, which gives 
enough time for nuclear decay photons to go through the 
SN ph otosphere. Compared to the models of Woosley et al. 
(19951, we notice that our model produces more 7 Be and 
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Figure 1. The upper panel represents the light curve data of 
SN1987A against our model for different time delays. The lower 
panel depicts the SN1987A data with the tdelay = 4 model. The 
SN1987A data comes from |Suntzeff et al.|fl99l| . The time delays 
represented in the upper panel are tdelay = 3, 4 and 9 days. In 
the upper panel we also depict the original M^i = 0.1 Mq target 
without the spallation, where luminosity contributions of 57 Co 
and 60 Co were added artificially, and their values, 5.5 X 10~ 3 Mq 
and 1.14 X 10 — 5 Mq, respectively, were taken from |Woosley et al.| 
(1995). The points are the SN 1987A data and the solid lines 
indicate our model. The y-axis is the logarithm of luminosity in 
units of erg/sec. In the lower plot, we also include the individual 
luminosity contributions of each isotope, where the contributions 
of 57 Co and 60 Co were artificially constructed in identical fashion 
to the upper panel. The masses of 56 Ni, 22 Na and 44 Ti, are actual, 
spallation products of our model. For tdciay = 4 days, the values 
are 0.99 M©, 5.50 X 10~ 5 M©, and 1.20 X 10~ 4 M©, respectively. 



22 Na. We also notice that 44 Ti, 56 Ni, and 26 Al are in the 



same order of magnitudes as other models (Woosley et al 
|1995[ ). Finally, it seems we can constrain production of Na 
and 6 A1 at higher tdeiay, where 22 Na production ends at 
tdciay ~ 7 and 26 A1 production at tdciay ~ 8 days. 



4 DISCUSSION AND PREDICTIONS 

• 44 Ti production as a dsQN specifi city and SN19 87A We 



argue that the rarity of observations ( The et al. 2006 1 of Ti 
might be related to the fact that 4 Ti is not produced by 



standard c c-SNe, but by dsQNe. In a recent paper (Ouyed 
et al.|201l[ ), we arg ued that CasA's 44 Ti lines were produced 
by a dsQN through a significant 56 Ni depletion. Although, 
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Figure 2. The panel represents the fractional abundance of 
atoms of a particular atomic mass number produced by our spal- 
lation model at a tdelay = 4 days, where the x-axis represents the 
mass number A, and the y-axis represents the fractional proba- 
bility of nucleus of mass number A. The panel depicts the yields 
for spallation of the inner ~ 1.0 Mq of the SN ejecta. 
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Figure 3. The panel represents the mass yields of 56 Ni, 44 Ti, 
22 Na, 26 Al and 7 Be, where x-axis represents tdelay and the y-axis 
represents the mass in units of Mq . The panel depicts the yields 
for spallation of the inner ~ 1.0 Mq of the SN ejecta. 



it is probably possible to avoid the excessive Ni depletion 
in our original model by adjusting some parameters, and 
yet still predict a 44 Ti mass yield that is compatible with 
observations, it seems more natural and elegant to extend 
our model with mixing. Through mixing, we can predict 
very minimal destruction of 56 Ni (Fig. I3J) because 56 Ni is 
spread throughout the whole envelope as opposed to the 
assumption where it is all concentrated in the inner layers. 
If 56 Ni is spread throughout a wider volume, then less 56 Ni 
atoms will be available in the inner layers as targets for 
QN neutrons. This minimal depletion leaves almost all the 
~ 0.1 M Q of 5B Ni intact. In Fig [j] we plotted the data of 
SN1987A against our model, where S N1987A is a SN where 
44 Ti has been measured i ndirectly (|Motizuki fc Kumagai 
2004J ISuntzeff et al.||1992| IWoosley & Hoffman||1991|). At 



lower time delays, QN energy is mostly spent in PdV work 
| Leahy fc Ouyed|2008 l, therefore at a time delay of 4 days, 
our model could reproduce a similar luminosity to SN1987A 
. Current observations of SN1987A point to thorough mixing 



of SN layers ( Mueller et al. 1991 1, which compares well to the 
assumptions about mixing done in our model. Recent papers 
have provided alternative arguments for the possibility of 



SN1987A compact remnant being a Quark Star (Chan et al 
[20091 . 

• 'Be abundance In Fig. [| and § a very large abun- 
dance of 7 Be is produced compared to current cc-SN models 
| |Diehl fc Timmes||1998| jWoosley et al.||1995[ ). Current core 
collapse SN models don't produce enough 'Be to reach the 
photosphere before 7 Be decays. However, our model predicts 
~ 10~ 3 - 10~ 2 M for 7 Be, which could be enough so that a 
significant amount of it reaches the photosphere and there- 
fore be detected by 7-ray telescopes. The 487 keV photon 
released by the /3-decay of 7 Be into 7 Li, is in the detectable 
energy range of current 7-ray telescopes (INTEGRAL, NuS- 
TAR, etc.). The large abundance of 7 Be is due to the spal- 
lation of mixed O and C in the inner SN layers. This high 
abundance in our model indicates that only core collapse SN 
that evolve into dsQN will have a detectable gamma signa- 
ture of 7 Be. Currently, Novae are the only possible sources 
of detectable, 7 Be photons because standard cc-SN are too 
optically thick ( |Diehl fc Timmes||1998| ). However our dsQN 
model can produce an alternative point source of 7 Be pho- 
tons. Finally, such a massive amount of 7 Be should create a 
signature in the bolometric light curve. 

• 22 Na abundance A small 22 Na contribution is predicted 



in late time core collapse SN light curves (Woosley et al. 
|1995[ ). However, our spallation model generates a much 
larger amount of 22 Na than expected. This large yields of 
22 Na are due to the spallation of 56 Ni into lighter nuclei. 
Current models ( jWoosley et al.|l995| > give a ~ 1O~ 6 M for 
22 Na, while our spallation model can produce as much as 
~ 10" 4 M© (Fig||J. Furthermore, 22 Na has a half-life of 2.6 
yr, and decays into an excited state of 22 Ne, releasing ener- 
getic radiation (Diehl & Timmes 1998). That large amount 
generates enough radiation to create a noticeable signal in 
the late time light curve, as pictured in Fig[l] The decay of 
22 Na emits a 1.275 MeV 7 line that can be detected by cur- 
rent 7-ray telescopes. These signals of 22 Na abundant dsQN 
cannot be reproduced in current models for cc-SNe. 

• Mixing Our model makes the simplification that 0,C 
and 56 Ni are in a perfect, homogenous mix. Our simplifica- 
tion allows us to model the spallation targets with statisti- 
cal weights derived from the isotopic abundances. However, 
Rayleigh Taylor instabilities induce mixing by fast moving, 
mushroom-like structures, structures that are lost when as- 
suming a perfect, homogenous mixture. Furthermore, stud- 
ies of CasA's remnants point to considerable asymmetry in 



the detonation (Fesen et al. 2006 1, which seems to contradict 



the assumption of perfect homogenous mixing. However, the 
main ideas of our model, that 44 Ti can be produced from 
a minimal destruction of 56 Ni, that spallation of 56 Ni will 
lead to enrichment of 22 Na, and that spallation of O, and 
C will lead to an excess of 7 Be, seem sound. While other 
limitations o f our model exist, which are discussed in | Ouyed] 
|et al.| ( |201 1\ , we offer some sound predictions for 22 Na and 
'Be enrichment . 
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